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Abstract 

Surface structure of a restricted ballistic deposition(RBD) model is examined 
on a one-dimensional staircase with free boundary conditions. In this model, 
particles can be deposited only at the steps of the staircase. We set up recur- 
rence relations for the surface fluctuation width W using generating function 
method. Steady-state solutions are obtained exactly given system size L. In 
the infinite-size limit, W diverges as L" with the scaling exponent a = ^. 
The dynamic exponent /? {W ~ t^) is also found to be ^ by solving the recur- 
rence relations numerically. This model can be viewed as a simple variant of 
the model which belongs to the Kardar-Parisi-Zhang (KPZ) universality class 
{ctRPZ = ^, I3kpz = 5)- Comparing its deposition time scale with that of the 
single-step model, we argue that (3 must be the same as (3kpz/{^ — (^kpz), 

which is consistent with our finding. 

PACS numbers: 05.40.-hj, 68.35.Fx, 68.55.Jk 
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I. INTRODUCTION 



Recently, dynamics of growing surfaces or interfaces has been one of the most interest- 
ing problems in surface science. Crystal growth, dielectric breakdown, fluid displacement in 
porous media, vapor deposition, spray painting and coating, biological growth and electrode- 
position are only a few examples which are important from both theoretical and practical 
point of view 0J^. Stochastically growing surface exhibits nontrivial scaling behavior and 
evolves to a steady state without any characteristic time or spatial scale. This observation 
has led to the development of general scaling approach for describing growing surface which 
exhibits a self-affine fractal geometry [^. In particular, the dynamic scaling approach [^,|^ 
has been applied to the study of a variety of theoretical models of growing surfaces. 

The important feature of growing surface is the roughness of surface induced by stochastic 
noises and its roughness can be characterized by few scaling exponents. A fundamental 
question is how to classify the rough surfaces given dynamic rules of deposition processes. 
The simplest model is the random deposition model which has no interactions between 
neighboring columns. This model can be solved exactly for any system size. It produces 
a rough surface in which the surface width grows with the square root of time but never 
reaches a steady state for finite systems. When interactions are introduced between columns, 
the surface width for finite systems saturates in the long-time limit. The evolution of the 
surface width W is described by a dynamic scaling relation 



where L is the linear dimension of the system, t is the time, and the scaling function / has 
the asymptotic behaviors f{x) ~ x^, with (3 = a/z, for x <^ 1, and f{x) ~ constant for 
X ^ 1. Consequently the limiting behavior of the width are 



Many stochastic models with intercolumnar interactions have been introduced and ex- 
tensively studied. One simple example is the random deposition model with surface diffusion 
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0. A deposited particle diffuses on the surface and is incorporated into the deposit upon 
reaching a local minimum. This rearrangement process reduces the surface tension of the 
deposit and the surface becomes relatively smooth with scaling exponents a — | and P — \ 
in (1 + 1) dimensions. These values are obtained by numerical simulations and the exact 
solutions are not available as yet. It has been argued that this model can be described by a 
continuum equation with a noise source 

_ = uW'h + V, (3) 

where h denotes the deviation of the local height of the deposit from its average value 
and ?7 is a Gaussian random noise. The solution of eq.(|^) yields a = 1/2 and j3 = 1/4 
in (1 + 1) dimensions [0, which is consistent with numerical results. Many other models 
with different rearrangement rules have been studied by numerical simulations and also by 
solving corresponding continuum equations. However, both methods have their own weak 
points. Numerical simulations often do not reach the scaling regime within a reasonable 
amount of time. A continuum equation corresponding to a given stochastic model may be 
reasonably guessed but it is very difficult to prove that it correctly describes the stochastic 
model. For example, there has been a great deal of controversies over which continuum 
equation describes the Wolf- Villain model [If. It is important to find exact solutions for 
stochastic models but it is usually formidable. 

Some other stochastic models like the ballistic deposition model f^J^, the single-step 
model P, p^ , p!T| , and the restricted solid-on-solid model [0 behave differently from the 



particle-rearrangement models. Numerical simulations on these models provide the scaling 
exponents a — | and /3 ~ | in (1 + 1) dimensions. In these models, the lateral current of 
particles along the surface is not conserved, so evolution of the surface involves the lateral 
growth. It has been argued that all these models can be described by the Kardar-Parisi- 



Zhang (KPZ) equation [jl3 



^ = uV'h + ^{Vhr + V (4) 
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where the nonhnear term represents the lateral growth. Dynamic renormalization group 
calculations for eq.(^ yield a = 1/2 and /3 = 1/3 in (1 + 1) dimensions (KPZ universality 



class), which is consistent with numerical results. However, higher- dimensional cases are not 
settled down yet. 

In this paper, we present a nontrivial stochastic model which can be solved exactly for 
finite system sizes in (1 + 1) dimensions. This model is similar to the ballistic deposition 
model, but height differences between neighboring columns are restricted to be positive or 
zero: hij^i — hi > where hi is the column height at site i. The resulting surface should 
look like a one-dimensional ascending staircase with variable widths of terraces and variable 
heights of steps (see Fig.l). Particles can be deposited only at the steps of the staircase, 
otherwise dynamic processes are rejected. From the master equation of this model with 
open boundary conditions, we set up recurrence relations for the surface fluctuation width W 
using generating function method. In the steady-state limit, we are able to extract the exact 
solutions for system size L from which we obtain a = | in the infinite-size limit. We also find 
/? = I by solving the recurrence relations numerically. No statistical errors are involved in 
this estimate because the ensemble average has already been done in the recurrence relations. 
It is somewhat surprising at first glance because this model can be viewed as a simple variant 
of the single-step model (KPZ universality class) |5|JTT[1 by rotating the surface clockwise by 
45° but the values of the exponents are different from those of the single-step model. We 
argue that the dynamic exponent (3 must be the same as Pkpz/{^ — Pkpz) by comparing 
its deposition time scale with that of the single-step model, which is consistent with our 
finding. Considering this argument other way around, one can say that our results for these 
scaling exponent are the best estimates for the stochastic models in the KPZ universality 
class. 

The outline of this paper is as follows. In section II, we describe our model and discuss 
its relation to a particle dynamics model, a spin-exchange model and the single-step model. 
In section III, we define the generating function and set up the recurrence relations. In 
section IV, we find the exact solutions in the steady-state limit and present the numerical 



estimate for (3. Comparison with the single-step model is given in section V and we conclude 
in section VI with a brief summary. 

II. RESTRICTED BALLISTIC DEPOSITION MODEL 

Consider the ballistic deposition (BD) model ||^ in (1 + 1) dimensions. In the BD model, a 
particle falls down along a straight line and sticks to the surface of the deposit. We introduce 
a restricted ballistic deposition model (RBD) by allowing only deposition processes which 
preserve the morphology of an ascending staircase of the surface, i.e. 

ni{t) = hi+iit) - hiit) > for all z, (5) 

where hi{t) is the column height of site i at time t and rij(t) is the height of the step between 
i-th and {i + l)-th sites. 

We start with a surface of horizontal length L with unit step heights at every possible 
steps, i.e. ni{0) = 1 for z = 1, 2, ■ ■ ■ , L — 1, equivalently hi{0) = i for i = 1, ■ ■ ■ , L. (The 
reason why we take this particular configuration as an initial state will be discussed in section 
IV). First, select a site i randomly and drop a particle along a straight line. The particle 
becomes part of the deposit when it contacts a particle in the deposit, only if the resulting 
step heights are all positive or zero (eq.(|^)). Otherwise, the particle is rejected(see Fig.l). 
When a site 1, L) is chosen at time t, the deposition process can occur for riiit) > and 
the rejection process for nj(t) = 0. Then and are updated as 

ni{t + At) = 0, 

ni_i(t + At) = ni_i{t) + ni{t), (6) 

where At is the time elapsed during one deposition attempt. Notice that eq.(|) is valid for 
both deposition and rejection processes. This is the key point of obtaining rather simple 
recurrence relations to be solved exactly. 

We use free boundary conditions. When a site i = 1 is chosen, the deposition and 
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rejection process can be described by the first equation of eq.(|D. At site L, only the 
deposition process can occur (no rejection) such that 

nL-i{t + At)=nL-i + l. (7) 

This equation can be identified as the second equation of eq.(^) by defining = 1. (This 
is useful to write simple recurrence relations in section III). 

The RBD model can be mapped on a one-dimensional particle dynamics model with 
directed diffusion and mass-conserving coalescence processes. Dynamic variable Ui represents 
the mass of the particle at site i. Eq.@ implies that a particle of mass nj(t) at site i diffuses 
to the left and stays there if nj_i(t) = (no particle at site i — 1 at time t) or coalesces with 
a particle of mass nj_i(t) at site i — 1 if nj_i(t) > 0. At the boundaries, there are input 
(source) and output (sink) of particles. Particles of unit mass come into the system from 
the right side of site L — 1 and particles of arbitrary mass leave the system to the left side 
of site 1 as shown in Fig.l. 

One can map this RBD model onto a variant of spin-exchange models, in the same way 
as the Toom model P^JT^ in the low-noise limit is mapped onto one of the spin-exchange 



models. Details of this mapping can be found in reference (16). Our model does not belong 
to general Ai^'^^ spin-exchange models(see Appendix A of |16|). This aspect will be discussed 
further elsewhere. 

One can find close resemblance between the RBD model and the single-step model if not 
sticky particles but smooth particles are used in the RBD model. Now a particle falls down 
along a straight line and becomes part of the deposit when it hits the ground. Of course, the 
surface of the RBD model is rotated by 45° from the surface of the single-step model. Both 
models restrict depositions only at the steps (or valleys on the 45°-rotated surface, see Fig.2). 
Unfortunately our model with smooth particles does not yield simple recurrence relations in 
contrast to the RBD model with sticky particles. Dynamics with smooth particles can not 
be represented by simple equations like eq.(|^) but equations with conditional statements. 
However, we argue in section V that these two models exhibit essentially the same scaling 
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behavior if time-scale difference is taken into account. 



III. GENERATING FUNCTION AND RECURRENCE RELATIONS 

In this section, we write the master equation for the evolution of the surface in the RBD 
model with general inhomogeneous deposition rates. The deposition rates are parameterized 
by K2, ■ ■ ■ 5 ^l} where Ki is a relative deposition rate at site i with normalization = 
1. The RBD model with ordinary homogeneous deposition rates is restored when all /tj's 
are replaced by A surface configuration can be characterized by a state vector n = 

(ni, 77-2, ■■■ , "^L-i) if we are not concerned about the overall average height but are interested 
in relative heights only. We define -Pt(n) as a probability to find a configuration n at time 
t. Then the master equation can be written as 

oo 

-Pt+At(n) = Ki^ni.o Pt{m, n2, ■ ■ ■) + KLO{nL-i)Pt{ni, ■ ■ ■ , ul-i - 1) 

m=l 

L-1 rii-i 

+ J2 f^fiiiT"i-i)K,o J2 ^ti- ■ ■ , ni_2, Ui^i - m, m, rij+i, ■ ■ ■) 

i=2 m=l 
L-1 

+ J2 t^iSn^fiPtin), (8) 
1=1 

where 5 is a Kronecker delta function and 6{x) = 1 for x > 1 and zero otherwise. The first 
three terms in the right hand side of the above equation describe the deposition processes at 
the boundaries and inside the boundaries. The last term describes the rejection processes. 
Consider the generating function defined as 

Gtizu ^2, • • • , ^L-i) = E ^I'^T ■ ■ ■ 4--i'Pt{n). (9) 

ni>0 

Using eq.@, one can obtain the generating function at time t + At in terms of Gt 

Gt+At{zi, ■ ■ ■ , Zl-i) = HlGt{l, Z2, ■■■) + KLZL-lGt{zi, ■■, Zl-i) 

L-1 

) Zi—ly Zi—iy Zi^i, ' ' ')• 

(10) 

i=2 

The ensemble averages of step heights (particle masses) and two-point (mass- mass) correla- 
tion functions are obtained by differentiating the generating function and setting all Zj's to 
be 1 as 
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{ni)t = -K— 

OZi {z}=l 

d ( dGi 



dzi V dzj 



(11) 

{z}=l 



Thus the dynamics of the step heights (masses) and two-points correlation functions are 
given as 



{ni)t+M = 


(1 


- i^i){ni)t + 




{ni)t+At = 


(1 


- '^i){ni)t + 2Ki+i{nini+i)t + Ki+i(n-+i)t, 




n'ini+i)t+At = 


(1 


- Ki- Ki+i){nini+i)t + Ki+2{n,^i+2)t, 


(12) 


{ninj)t+At = 


(1 


- Ki- Kj){ninj)t + Ki+i{nj+inj)t + Kj+i{ninj+i)t. 





where i, j — 1, ■ ■ ■ , L — 1 with |i — j | > 2 and ul is set to be 1. 

First, we consider the steady-state only. The ensemble-averaged masses in the steady 
state, (nj)oo, satisfy the recurrence relations 

f^Mijoo = <x> = • • • = Kl(^^l)oo = I^L, (13) 

which leads to 

(n^oo = -. (14) 

In the homogeneous case, (ni)oo = 1 at any site i, so the ensemble-averaged surface looks 
like an ascending regular staircase with unit step heights and unit terrace widths in the 
steady-state (long-time) limit. 

We also find recurrence relations for two-point correlation functions in the steady state 

as 

QiL = Qi+1,L, 
Qii — '^iQi,i+l + Qi+l,i+l, 
Qi,i+1 = 1fi,i+lQi,i+2) (15) 

Qij — UijQi^ij -\- VijQij^i, 



where Qi-,-, Uij, Vij, and Wi are defined as 

Qij = Ki{ninj)^, (16) 

Uij = ^ , , Vij = Uij^^^, and Wi = 2^^, (17) 

with Qji = ^Qij- It is easy to solve the above recurrence relations numerically for finite 
values of L. Start from the known value oi Qll = i^l- Using the recurrence relations 
repeatedly, one can find the values of Ql~i,l, then Ql-2,l and Ql-i.l-i, then Ql-3,l and 
Ql-2,l~i, then so on. Qij is determined by Qi'ji's only for i' + j' > i + j. This brings out 
an idea of mapping these recurrence relations to a directed random walk problem. 

Consider an upper diagonal half of a L x L square lattice (Fig.3). We assign appropriate 
weights; Uij, Vij, Wi, 1, or on the bonds connecting neighboring sites. Now imagine that a 
walker starts from a lattice site (L, L) and moves only left or down (directed random walk). 
There are numerous paths along which the walker can reach a site A typical path is 

drawn in Fig.3. Paths going through a lower diagonal half of the square lattice are excluded. 
The weight to a path is given as the product of bond weights along the path. Then, Qij can 
be calculated by summing up path weights for all possible paths. So the formal solutions 
for Qij are given as 



Qij = Qll Yj Yl {uij, Vij, Wi, 1,0). (18) 

all a path 

possible paths 

This is similar to the first passage problem when all bond weights are set to be equal. We 



are able to extract exact solutions in the case of homogeneous deposition rates (i.e. Ki = 1/L 
for all site i). 

IV. EXACT SOLUTIONS AND CRITICAL EXPONENTS 

We consider the homogeneous case only. Then Uij = u = 1/2, Vij = v = 1/2, and 
Wi = w = 2. And Qij = Qji. For convenience, we define normalized two-point correlation 
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functions Qij = Qij/QiL where Qll = 1/L. From eg. ([T8|) , Qij is the sum of path weights 
over all possible directed paths from (L, L) to which do not go through a lower diagonal 
half of the L x L square lattice. Some of these paths have the same path weights but some 
do not. So it is necessary to distinguish the paths with different path weights. 

First, consider the case i < j {j < L — 1). Paths going through a bond connecting the 
points (£, L) and {i,L — 1) > i) must have the same path weights, 

• I \ L+e-i-j 



w(^,j;£)=i^-v-v-*= 



(19) 



The number of such paths, M{i,j\€), is equivalent to the number of directed paths from 
(£, L — 1) and to which do not touch the diagonal line. Paths touching the diagonal 



line do not contribute to Qij for i < j. Using the reflection principle [0, the number of 
paths can be easily obtained as 



/ 



- 1 - ? ^ ( i- 



V 



V 



t+L-l-j 
L-l-i 



L-l-j 

where (■ ■ ■) is a combinatorial factor. Introducing new site variables for convenience {p 
L — i — 1, q = L— j — l,r = L — £ — 1), the normalized correlation functions are 



(20) 



L-2 



(21) 



1 ^ 

^ r=l 



p + q-r 



V 



-r ]^ q 

' 2 ^ 



r=l 



p 



J 



\ \ P+1- 
2) 



This equation can be rewritten as 



Qij - 2 

where k = p — q = j — i>0 and 



q p 
m=l—k m=l+k 



(22) 



Z2nim] 



^ 2n — 



n 



2n—m 



(23) 
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which can be interpreted as a probabihty that a random walker in one dimension returns to 
the starting point exactly m times up to 2n steps |jT^. Using identities (eqs.(Al) and (A2) 

in appendix), eq.(|2^) simplifies to 

fc-i p 
Qi, = J2 Z2p{m) = 1 - E ^2p(m). (24) 

m=0 m=k 

The diagonal elements Qa can be written, using eg. (p!5D , as 

L-l 

Qii = 1 + 2 E Qjj+i 
j=i 

= l + 2(p+l)Z2(p+i)(0), (25) 

where eg. (|2^) and the identity of eg. (A3) are utilized to derive the second eguality. Egs.(piD 
and (^) form a complete set of exact solutions for mass-mass correlation functions in the 
steady state. 

Now we are ready to calculate the fiuctuation properties of the surface in the RBD model. 
The surface height at site i, hi{t), can be written in terms of step-height variables %(t) as 

kit) ^Y.^ At), (26) 

where i = 2, ■ ■ ■ , L and the reference height hi{t) is set to be zero. We consider two important 
fiuctuations of the surface; the step-height (mass) fiuctuations and the height fiuctuations. 

Mfmt)^{n^)t-{n,)l 

W^{L,t)^{hl,),-{h,^,f, (27) 

i 

= E [{i^jr^j')t - {nj)t{nj')t] ■ 

Notice that the fiuctuations are site-dependent due to the lack of the translational symmetry 
of our RBD model. So it is useful to consider the fiuctuations averaged over all sites. The 
average mass fiuctuation in the steady state is 

M\L, oo) = -^J2M^ = hz^LiO), (28) 
^ - ^ i=i -3 

where eg. (|25|) and the identity of eg.(A4) is used for the second eguality. For large L, 
grows like L^"' with exponent a' = 1/4. It is a little tricky to calculate the average height 
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fluctuation W"^ in the steady state. However, using few identities of eqs.(A4) — (A7) and 
eqs.(^) and (|25|), we find that it is simply related to the average mass fluctuation as 

W\L, oo) = 2L- M\L,oo) . (29) 

For large L, W"^ grows hke L^" with exponent a = 1/2. 

Time-dependence of the fluctuations, and W"^, can be investigated numerically, using 
the time-dependent recurrence relations in eq.(|l2]). We take a surface with unit step heights 
at every possible step (72^(0) = 1 for all i) as an initial configuration which is equivalent to the 
ensemble-averaged surface configuration in the steady-state limit ((?T.j)oo = !)• If some other 
surface configurations like a flat surface are taken as an initial configuration, the dynamics 
of the growing surface involves two different time scales; one time scale associated with 
the evolution of the ensemble-averaged surface into the steady-state configuration and the 
other time scale associated with the development of surface fluctuations. The former time 
scale is trivial and uninteresting. Moreover, the scaling exponent (3 in eq.(|l|) characterizes 
the surface-fluctuation time scale. We can single out the surface-fluctuation time scale by 
choosing the above initial configuration. Starting with (nj)o = 1 and (njnj)o = 1, we solve 
the time-dependent recurrence relations, eq.([T2|), iteratively. In this paper, we consider the 
homogeneous case only; Ki = 1/L. The results for the average mass fluctuation and the 
average height fluctuation W"^ versus time t are plotted in Fig. 4 (log-log scale) for system 
size L = 400. We find that ~ t^/?' ^j^j^ ^/ ^ 0.253(5) ~ 1/4 and ~ t^/? ^j^^ 
P = 0.500(1) ~ 1/2. In both cases, the dynamic scaling exponents, z = a/ (3 and z' = a'/P', 
are found to be 1 {z = z' = 1). 

V. TIME SCALES IN THE RBD MODEL AND THE SINGLE-STEP MODEL 

The RBD model can be viewed as a simple variant of the single-step model by rotating 
the surface by 45°. As discussed in section II, the RBD model with smooth particles instead 
of sticky particles is equivalent to the single-step model except for boundary conditions. 
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Depositions occur only at the valleys of the 45°-rotated surface. With smooth particles (or 
in the single-step model), the height of the valley increases by one unit at the deposition of 
a particle. However, with sticky particles in the RED model, the heights of all sites between 
the valley and the hill increase by one unit (see Fig.2). In average, the deposition of a 
particle in the RED model is equivalent to the depositions of many particles between the 
valley and the hill in the single-step model. We argue that the number of depositions in the 
single-step model corresponding to a single deposition in the RED model is proportional 
to the average surface fluctuation width W. In order to get the same morphology of the 
surface in average, one should wait for much longer time ts in the single-step model than in 
the RED model. The elapsed time tg in the single-step model is related to the elapsed time 
t in the RED model as 

ts ~ Wt. (30) 

The surface fluctuation width of the single-step model Wg at time tg is equivalent to that of 
the RED model at time t in the infinite-size limit; 

W{oo,t) = Ws{^,ts). (31) 

The single-step model belongs to the KPZ universality class with scaling exponents axpz = 
1/2 and Pkpz = 1/3. Using Ws{(x>,ts) ~ tf^^^, eq. (pip becomes 

iy(cx),t) ~ (m)^^^^, (32) 

therefore 

W{oo,t)^t^, (33) 

with 

J- — Pkpz ^ 

It is shown that the scaling exponent /? is 1/2 for the RED model by simply comparing 
the deposition time scales in the RED model and the single-step model, which is consistent 
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with our result in the previous session. It imphes that these two models exhibit essentially 
the same scaling behavior if the time-scale difference is taken into account. Therefore the 
precise measurement of the exponent j3 in the RBD model (which is done in section IV) 
provides the good estimate for the scaling exponents of the stochastic models in the KPZ 
universality class. 

We perform Monte Carlo simulations to confirm our argument about the time-scale 
difference. We take the initial configuration as a surface with ^^(O) = 1 for all i with 
horizontal length L = 500. After one deposition attempt on the average per lattice site 
(one Monte Carlo Step) , the time is incremented by one unit in the time scale of the RBD 
model. In order to use the time scale of the single-step model ts, the time increment for an 
actual deposition at site i must be multiplied by the step height ni{ts). In this time scale, 
we run simulations until tg = 5, 000 and the number of independent runs is 3,000. The 
double-logarithmic plot (Fig. 5) for the surface fluctuation width W versus time tg shows 
a straight hne with the slope Ps — 0.338(6), which agrees with the KPZ value 1/3. This 
result confirms our previous argument about the time-scale difference between the RBD 
model and the single-step model. Moreover, the mass fluctuation M in this time scale also 
shows the scaling behavior with exponent P'^ — 0.168(4) ~ 1/6. In both cases, the dynamic 
exponents Zs and z'g take the KPZ value 3/2. This result strongly supports our suggestion 
that the RBD model exhibits essentially the same scahng behavior as the models in the KPZ 
universality class if the proper time scale is employed. 

VI. SUMMARY AND CONCLUSION 

We have studied a restricted ballistic deposition (RBD) model on a one-dimensional 
staircase with free boundary conditions. Using generating function method, the exact solu- 
tions for the surface fluctuation and the mass fluctuation were obtained for any system size 
in the steady-state limit. The RBD model is one of a few nontrivial stochastic models which 
can be solved exactly. We also examined time-dependent solutions by solving the recurrence 
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relations of two-point correlation functions numerically. Prom these solutions, we extracted 
the values of the scaling exponents a and /3, both of which turned out to be 1/2. 

The RBD model is quite different from the ordinary ballistic deposition (BD) model 
because the RBD model does not allow the down steps. Indeed, the scahng behavior of the 
RBD model differs from the BD model which belongs to the KPZ universality class where 
OLKPZ = 1/2 and Pkpz = 1/3. However, the RBD model can be viewed as a simple variant 
of the single-step model which also belongs to the KPZ universality class. These two models 
exhibit a very similar surface morphology. And their dynamical processes are not much 
different in average except for the difference between deposition time scales. 

We suggested that the major difference in the scaling behavior of the RBD model and the 
single-step model may disappear if the time-scale difference is properly taken into account. 
So the static exponent a should be equivalent in both models but the dynamic exponent /3 
may be different. We explored the relation between two exponents as P — Pkpz/{^— Pkpz)-i 
which is consistent with our results. Monte Carlo simulations confirm our argument about 
the time-scale difference. 

The generalization of the RBD model into higher dimensions may be quite interesting, 
partly because it may serve as a useful alternative model for the indirect investigation of the 
KPZ universality class in higher dimensions. Also our generating function approach allows us 
to investigate the models with inhomogeneous deposition rates. Effect of the inhomogeneity 
on the surface morphology and the scaling behavior in the RBD model is currently under 
study. 
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APPENDIX 



Some useful identities are listed in this appendix. 

n 

Z2n{m) = 1. 



m=0 



fe-1 k-l 

^ Z2n(m) = Z2(n-k){-m - 1). 
m=0 m=0 

5]Z2„(0) = 2pZ2p(0). 

n=0 

|^nZ2„(0) = ^p(p- 1)^2^(0). 

n=0 

I: n2Z2„(0) = - l)(3p - 1)^2^(0). 

n=0 -"-^ 

mZ2n(m) = (2n + 1)^2^(0) - 1. 

m=0 

n 

^ m'Z2„(m) = -3(2n + l)Z2n(0) + (2n + 3). 

m=0 
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Figure captions 

Fig.l : Dynamic processes of the RBD model. Depositions of particles can occur only at 
the steps of the surface. The corresponding particle dynamics model is also shown. 

Fig. 2 : Deposition processes of (a) the single-step model and (b) the RBD model on the 
45°-rotated surface. 

Fig. 3 : Recurrence relations are mapped to a directed random walk problem in an upper 
diagonal half of a L x L square lattice. The thick line represents a typical path along 
which a directed walker can start from {L,L) and end at 

Fig. 4 : Log- log plots of the average mass fluctuation M^(«) and the average height fluctu- 
ation W'^{o) as functions of time t. These data are obtained by solving the recurrence 
relations, eq.(]l2D. The slopes of the lines are 2(3' = 0.50 and 2(3 = 1.00. 

Fig. 5 : Log-log plots of the average mass fluctuation M^(») and the average height fluctua- 
tion W'^{o) as functions of time ts- These data are obtained by Monte Carlo simulations 
using the time scale of the single-step model. The slopes of the best-fitted hues are 
2P' = 0.336 and 2/? = 0.676. 
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